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STATIONARY MEASURE OF THE DRIVEN TWO-DIMENSIONAL 
g-WHITTAKER PARTICLE SYSTEM ON THE TORUS 

IVAN CORWIN AND FABIO LUCIO TONINELLI 


Abstract. We consider a g-deformed version of the uniform Gibbs measure on dimers on the 
periodized hexagonal lattice (equivalently, on interlacing particle configurations, if vertical dimers 
are seen as particles) and show that it is invariant under a certain irreversible q- Whittaker dynamic. 
Thereby we provide a new non-trivial example of driven interacting two-dimensional particle system, 
or of (2 -I- l)-dimensional stochastic growth model, with explicit stationary measure. We emphasize 
that this measure is far from being a product Bernoulli measure. These Gibbs measures and 
dynamics both arose earlier in the theory of Macdonald processes [7]. The q = 0 degeneration of 
the Gibbs measures reduce to the usual uniform dimer measures with given tilt [12] , the degeneration 
of the dynamics originate in the study of Schur processes um and the degeneration of the results 
contained herein were recently treated in m- 


1. Introduction 


Irreversible Markovian dynamics on two-dimensional dimers or interlaced particle configurations 
(see Figure are closely related to driven interacting particle systems as well as random surface 
growth models in (2 -|- l)-dimensions. It is a challenge to hnd local irreversible dynamics whose 
invariant measures (on the torus or on the inhnite lattice) are likewise local and explicit. Knowledge 
of a dynamic’s invariant measures can be useful in establishing its hydrodynamic / fluctuation 
theory and in understanding general properties (i.e. universality classes) of two-dimensional driven 
systems. 

The state space for the Gibbs measures / dynamics we consider is that of interlacing particles on 
a discrete torus. Each particle interlaces with two particles above it, two particle below it, and has 
two neighbors at the same row (see Figure [^. For a particle p we let Ap, Bp, Cp, Dp, Ep, Fp denote 
the absolute value of the horizontal distances between the particle and its neighbors starting with 
the right neighbor on the same row and going clockwise (actually, the values of Ap, Bp, Dp, Ep are 
this distance minus 1 - see Section 2.1 for precise definitions). In continuous time, each particle p 
jumps to the right by one lattice space according to an exponential clock of rate 


(l-g^^)(l-g^^+^) 

(1 - g'C'p+l) 


where q G [0,1). This dynamic is local and irreversible. Moreover, it preserves the interlacing of 
particles - particle p cannot jump past the particle below and to its right since the factor 1 — q^p = 0 
in that scenario; and if a particle p jumps so as to pass its upper-right neighbor, then that neighbor 
is immediately pushed right by one space since the denominator for that neighbor’s jump rate is 
zero, hence it jumps right at an inhnite rate. 

The main result of this paper. Theorem provides an explicit two-parameter family of transla¬ 
tion invariant Gibbs measures on the torus which are invariant for this dynamic. One parameter 
represents the number of particles per row, and the other is topological, related to how much the 
particles rotate right between rows (in terms of dimers, the parameters relate to the slopes of the 
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height function). The Gibbs measures are local and proportional to 


n 


( g ; q)ap 

{Q:Q)Bp{q;q)cp 


where the product is over all particles p and where {q]q)n = (1 ~ g)(l ~ g^) • • • (1 — q"') is the 
g-Pochhammer symbol. These dynamics are irreversible (i.e. driven); the proof of the invariance 
of the Gibbs measure involves unexpected cancelations of contributions from every particle. For 
q = 0, the dynamic reduces to the one introduced by A. Borodin and P. L. Ferrari [5] and the 
Gibbs measures reduce to uniform measures over the configurations of interlacing particles. Their 
stationarity on the torus was proved in |19) and the argument is simpler than in the q G (0,1) case. 

In general, there is no hope of computing explicitly the stationary measures of a driven particle 
system. In the lucky cases when it is possible, often the invariant measures turn out to be of 
product type. This is the case, for instance, for the ASEP on Z'^, d > 1 and certain modifications 
thereof m , zero range processes mass transport models [9], and certain discretizations of 

the one-dimensional KPZ equation including directed polymer models [miiTiiniE]. In contrast, 
our stationary measures vr are far from being product: for instance, for g = 0 it is known that they 
show power-law decaying spatial correlations and the same presumably holds for q ^ 0- 

The dynamics and Gibbs measures we study in this work were initially considered in the study 
of Macdonald processes (in fact, Whittaker processes) [3 Chapter 3], though our results do not 
rely at all on this technology. Still, let us briefly explain the origins. The g-Whittaker process 
is a measure on interlacing triangular arrays with 1 particle on the first (bottom) row, 2 on the 
second, up to N on the top (N^^) row. Given the configuration of particles on the top row, the 
g-Whittaker process enjoys the property that the measure on the remaining N — 1 lower rows is 
given by the above described Gibbs measure. The dynamics we described were also introduced in 
[3 Chapter 3] and play nice with such measures on triangular arrays. In particular, if one starts 
with an interlacing triangular array satisfying the Gibbs property, then after running the dynamics 
for a fixed time, the resulting triangular array also enjoys the same Gibbs property (of course, the 
top row will have moved to the right). This fact follows from results in [3 Section 2.3]. 

While our main result does not follow from the results of [ 3 , it was partially inspired by it. 
After all, if the Gibbs property is preserved on these triangular arrays, it would seem reasonable 
that it should also be preserved on the periodized version of the model. In [3 a more general 
class of dynamics and Gibbs measures are discussed which are inhomogeneous with row-dependent 
parameters Ofc > 0. Our approach applies equally well in that more general context and so we state 
and prove our main theorem for general (it reduces to the case described above when all = 1). 
There actually exist many types of continuous time Markov dynamics on triangular arrays which 
preserve the Gibbs property [8], however the other known dynamics do not easily admit periodic 
versions (i.e. they are not translation invariant). There is also a discrete time version of the above 
introduced dynamic described in [3 , however it is not as simple and we do not pursue studying its 
periodized version. 

The other inspiration for our result is the recent work |19) regarding the g = 0 case of this model 
(that work was inspired by the Schur process work of O [6] ). The g = 0 case of the dynamics is 
known to lie in the (2-|-l)-dimensional anisotropic Kardar-Parisi-Zhang universality class [5l[6lll9|. 
Based on the approach developed in [19], the results of the present paper can be seen as a step 
towards extending this universality to g G (0,1). In the g = 0 case, the next step in [T9| is to extend 
the invariance of Gibbs measures from the finite torus to infinite volume. A crucial ingredient used 
in this extension was the fact that, for g = 0, the infinite-volume Gibbs states are known and have 
an explicit determinantal structure and GFF-like height fluctuations [12] • All of this is missing in 
the g G (0,1) case, so at present, the extension of our result to the infinite lattice and the proof of 
its KPZ class behavior is an open problem. 
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Recently, Borodin-Bufetov [Ij considered another deformation of the Borodin-Ferrari particle 
system {q = 0 case of our dynamics) and they showed that the invariant measures are given by 6- 
vertex Gibbs measures. This deformation is different from the one we consider here (ours originates 
from g-Whittaker processes [7] and the other from vertex models [3j). Let us also mention that 
another example of (2 + l)-dimensional growth model with explicit (non-product form) stationary 
measure is the Gates-Westcott model [iniiii!. 

Outline. In Section]^ we introduce the state space for our Markov dynamics (in terms of dimers as 
well as interlacing particle configurations) and describe the class of Gibbs measures we will work 
with. Section introduces the periodized Whittaker dynamics. Section [^contains the proof that 
the Gibbs measures are invariant for the g-Whittaker dynamics. 


2. State space and Gibbs measure 

The dynamics we study can be described as a driven interacting particle system, or as dynamics 
of a dimer model on the (periodized) hexagonal lattice. While the former description may be more 
natural, the latter allows to get more directly certain statements, notably ergodicity of the Markov 
chain. We start by introducing the state spaces for our dynamics. 


2.1. Description as an interacting particle system. The particle process lives on the L x N 
discrete torus x 'Ll{N'L). The horizontal size is L and the vertical size is N. 

The particle configuration space will be denoted aiid depends on two integers 1 < 

mi < L and 1 < m 2 < N such that 


niilL -|- m2/N < 1 . 


( 2 . 1 ) 


At each site x = (xi,X 2 ) G there is at most one particle. On each row there are exactly mi 
particles. We exclude mi = 1 and mi = L to avoid trivialities. The parameter m 2 has a more 
topological nature and its meaning will be explained below. 

Particle positions are interlaced, in the following sense. The horizontal position of particle p is 
denoted Xp G Z/(LZ). Given any p (say on row i), we let pi,Pi denote its right/left neighbor on 
the same row (note that if mi = 2 then pi = p^). Then, we require that in row i — 1 there is exactly 
one particle, denoted p 2 , whose position satisfies 

^P2 ^ T T 2,..., Xpj } (^■2) 

and exactly one particle, denoted ps, satisfying 

Xp3 G {xp^ + Ij Xp^ -|- 2,..., Xp}. (2-3) 

See Figure Note that, automatically, in row i + 1 there are exactly one particle ps and one 
particle pq satisfying respectively 


‘'PS 


£ {Xp4,..., Xp 1 }, 


‘'P6 


G }xp ,..., Xp^ 1}. 


(2.4) 
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Figure 1 
allow Cp : 


The neighbors pi,... ,pq of particle p. Note that conditions (2.2), 
= 0 but they impose Bp + 1 := Xp^ — Xp> 1. 


(2.3) 
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We define non-negative integers Ap,... ,Fp as 


Ap — ^pi 

Dp — Xp 1, 


Bp — Xp2 Xp 1 , Cp — Xp Xp2 
Ep — Xp Xpg 1, Fp — Xpg Xp. 


(2.5) 


The particles pi,... ,p6 are the six neighbors of p, labeled clockwise starting from the one on the 
right. The definition of the dynamics will be such that the labels of the neighbors of a particle p 
do not change with time. 

Let ^L,N-,mi be the set of particle occupation functions, i.e. of functions p : i—)■ {0,1}, 

with mi particles per row, whose positions satisfy the constraints (2.2)-(2.4). The set 
decomposes into disjoint “sectors”: 




yyii ,m2 


( 2 . 6 ) 


as follows. Given any particle p, connect p to its up-right neighbor pg, then pg with its own up-right 
neighbor and repeat the operation until the path T thus obtained gets back to the starting particle 
p. Note that T forms a simple loop: otherwise, there would be a particle r which is reached along T 
from two different particles r', r". This is impossible, since both r' and r" would be the rg neighbor 
of r. Call G NU {0}, G NU {0} the vertical and horizontal winding numbers of T around 
the torus Tl,n- It is easy to see that N^, Ny are independent of the chosen initial particle p. It is 
possible to check (see Remark below) that 

Eh 

m2.= mi— (2.7) 

I\v 


is an integer and that it satisfies (2.1). The set is the subset of ^L,N-,mi such that m 2 = k. 

Each sector will remain invariant under our dynamics. 

Given q G [0,1) and a collection of real parameters ai,...,aN > 0, let vr be the probability 
measure on ^L,N-,mi,m 2 defined as 


’^(* 7 ) ■— ('^) ■ — 


ZL,N;m„n.2 ^ (q; g)« (9; ^ 


Cy 


{Q',q)a„ 


( 2 . 8 ) 


where {q; q)n = (1 — ( 7)(1 — q ^)... (1 — g”) and where r(p) denotes the label of the row in which p 
sits. 

Conditionally on the position of all particles except particle p, the law of the position of p is 
proportional to 


Cp Fp (g;g)Ap(g;g)F>p 

^r{pfr{pe) q)c^{q-, q)Eyiq; q)Fp ' 


(2.9) 


Note also that in (2.8) we could have for instance replaced {q; q)Bp{q] q)cp by (g; q)Ep{q] q)Fp and tt 
would be unchanged. Likewise, because Ylp-r{p)=k Cp + Bp = N — mi we could have replaced 

by ®p(p) any a and not changed f either. When g = 0 and all Ofc = 1, f reduces to the 

uniform distribution on ^L,N-,m^_,m 2 - 


Remark 1. Note that we are viewing ^L,N;mi,m 2 ® set of particle occupation functions and 
not as a set of positions of labeled particles; i.e. two particle configurations with the same particle 
occupation variable everywhere are considered to be the same. If instead we looked at configurations 
of labeled particles, every configuration in IlL,N-,mi,m 2 would actually correspond to mi different 
possible particle configurations, obtained by cyclically changing particle labels, simultaneously on 
every row. 
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2.2. Description as a dimer model. Before we introduce the dynamics, let us give the alternative 
description of the model in terms of perfect matchings (or dimer coverings) of the periodized 
hexagonal lattice Precisely, with reference to Figure]^ ^l,n has period L in the ei direction 

and period N in the 62 direction, and vertices are alternately colored black/white. 



Figure 2. The graph with L = 5, N = 4. The graph is periodic in both 
ei and 62 directions. Vertical dimers (in red) will be called “particles”. Here, 
v' = V + 2ei, v” = V + 62 + ei and v”' = v — 62 + ei. 


The number of dimers (i.e. of edges in the perfect matchings) equals LN, which is the number of 
white vertices. Given strictly positive integers ni,n 2 with ni + n 2 < NL, let ^L,N;ni,n 2 be the set 
of dimer coverings X with ni vertical dimers and n 2 north-west oriented dimers. It is well known 
that the number of vertical dimers is the same in each horizontal row, so that ki := ni/N G N, 
and similarly there are k 2 := n 2 /T G N north-west oriented dimers in each e 2 -oriented column. To 
avoid trivialities we will also assume that ki > 2 (recall that in the particle picture we required 
mi > 2 for the same reason). Another well-known fact is that the knowledge of which vertical 
edges are occupied uniquely determine the whole dimer configuration. A last observation is that 
the horizontal (i.e. ei) coordinates of vertical dimers are interlaced: given two vertical edges v and 
v' = V + kei, k > 0 such that there is a dimer at v, v' and no dimer at u -|- jei, j = 1,..., A; — 1, 
then there is exactly one value 0 < r < A: and one value 0 < s < k such that there is a dimer at 
v” = V + 62 + r 6 i and at v'" = v — 62 + S 6 i. See Figure 


Remark 2. If V 6 rtical dimers are called “particles”, then the correspondence between the particle 
picture of Section and the dimer picture should now be clear, see also Figure [^‘ for each particle 
configuration in IlL,N-,mx,m 2 there is a perfect matching of with certain values (A:i,A: 2 ), and 

vice-versa. It is obvious that mi = ki. It remains to be shown that also m 2 = k 2 , thereby proving 


also ( 2 . 1 ), since Lk 2 + Nki = ni+n 2 < NL by assumption. Given a nearest-neighbor path Cf^f/ 
from a face f to a face f and a perfect matching X, define 


HxiCf- 


■f) 


edesA: 


( 2 . 10 ) 


esc 




where the sum runs over the edges crossed by the path, Cg equals -1-1/ — 1 if e is crossed with the 

white vertex on the right/left and leex is the indicator function that e is occupied by a dimer in 

X. It is well known m that if f = f then, given X, Hx{Cdepends only on the horizontal 

and vertical winding numbers of Cf^f around the torus. Choose a path that moves only in 

the -\-ei direction. Then, we see that 

Hxicfif) = +ki. 


(2.11) 
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(2) 

If instead it moves only in the +62 direction, then 


Hx{C 


(2) N _ 


f^f> 


= -k2. 


( 2 . 12 ) 


(r) 

On the other hand, choose as follows: it starts from a face f just to the right of a vertical 

dimer p, it moves in direction +62 if this involves crossing no dimer, and in the direction +ei 
otherwise; the paths stops when it gets back to f. Let us say that a vertical dimer r “is visited by 
visits the hexagonal face just to the right of r. Note that the first visited vertical 
dimer is p, the second one is its up-right neighbor pQ, then the up-right neighbor ofpQ, and so on. 
These are just the particles visited by the path F built in Section \2.1\ In particular, we have seen 
that r forms a simple loop, so indeed come back to f. Also, has the same winding 

numbers Nh,Ny as F. By construction, we see that = 0 since the path never crosses 

any dimer. On the other hand, by (2.11 )-(2.12) we have 


0 = Hxicf_lf) = Nhki - Nyk2, 


(2.13) 


i.e. k 2 = m 2 , see (2.7). 



Figure 3. The correspondence between (a portion of) perfect matching and in¬ 
terlacing particle configuration. The mapping is made more evident by an affine 
deformation of the hexagonal lattice such that the axes ei, 62 become orthogonal 


In analogy with Section 2.1, for each vertical dimer p, let pi,...pQ be the labels of its six 
neighboring vertical dimers labeled in clockwise order starting from the one on the right. The non¬ 
negative (possibly zero) integers Ap, Bp, Cp, Dp, Ep, Fp defined in (2.5) are given in the coordinates 
of the hexagonal lattice as follows: if p is at vertical edge v, then 


Pi is at edge v + {Ap -\- l)ei; 

P 3 is at edge v — 62 — CpCi ; 

P 5 is at edge v -\- 62 - {Ep l)ei; 


P 2 is at edge u - 62 + {Bp -\- l)ei 
P 4 is at edge v — {Dp + l)ei 
P 6 is at edge u -|- 62 + FpCi. 
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3. PeRIODIZED g-WHITTAKER DYNAMICS 


We saw that the descriptions in terms of interlaced particles on the torus Tl,n or of perfect 
matchings of are equivalent, if we set {mi, m 2 ) = {ki,k 2 )- We will interchangeably use the 
former or the latter representation. 


Definition 1. Given a configuration rj £ draw a directed upward edge from any particle 

r to its up-right neighbor rg if Fr = 0 and a downward edge from r to r^ if Cr = 0 (in both cases, we 
draw an edge between particles in neighboring lines, with the same horizontal position). For each 
particle p let (resp. V~) be the set that includes p plus the particles that can be reached from p 
by following upward (resp. downward) oriented edges. 

Remark 3. From the assumption ni,n 2 > 0 it follows that \ V^\ < A^ —1 for every p. In fact, if say 
11^"^ I > N then the upward edges starting from p form a loop and actually \ Vj)~\ = N (in each row 
there is at most one particle p' with the same horizontal position Xpi as p). In this case, the path 
r defined just after ( |2.6[ ) has zero horizontal winding number, so from ( |2.7[ ) we have m 2 = n2 = 0. 
Therefore, the up/down arrows do not form loops and we can identify r(fi, the highest particle in 
, and r~, the lowest particle in V~. 


The dynamics we consider is a continuous-time Markov chain on ■ The updates consist 

in shifting by -\-ei all particles in one of the families . Such move happens with rate 


(1 - q^p){l - 

(1 _ qCp+i^ 


(3.1) 


where q £ [0,1) and ai,..., on > 0 are the same real parameters as in the definition of the Gibbs 
measure (2.8). Recall that r{p) is the row associated to particle p. Note that the rate is zero 
if Bp = 0. This prevents particles from overlapping after the move, see Figure]^ (b). Note also 
that after the move, the configuration is still in IlL,N-,ni,n2- Indeed, if = 1, shifting p by -|-ei 
corresponds to rotating three dimers around a hexagonal face (see Figure (a)), so ni,n 2 remain 
constant. In the general case, if Rp > 0 the shift of Vjfi by ei can be obtained by shifting -|-ei 
the particles of V^fi one by one, starting from the top one, i.e. r+. Also in this case, ni,n 2 are 
unchanged. 



Figure 4. (a): If |lp| = 1 then necessarily edge v is occupied and if Bp 0 then v' 
is also occupied. The allowed shift of ei then corresponds to the rotation of these 
three dimers, (b): A particle p with the associated set V.^, whose highest particle is 
r+. The shift of V.^ by ei should be forbidden, otherwise particles r+ and p' would 
overlap. Indeed, by the particle interlacing condition note that edges vi,V 2 ,vz are 
necessarily occupied by particles. In particular, Rp = 0 because U 3 is occupied, so 
that the rate (3.1) is zero. 
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Remark 4. For q = 0 and = 1, the dynamics was studied in [3 1 EH] (in this case, 
should be interpreted as 1 when Bp = 0 ). As we already remarked, in this case vr = TTL,N;mi,m 2 is 
the uniform measure over IE]; “particles” were associated to north-west oriented 

dimers, rather than to vertical ones as in the present work. With the convention of |19) . updates 
consist in a single particle jumping a distance n> 1 in the —62 direction, instead o/n > 1 particles 
jumping a distance 1 in the +ei direction. 

A preliminary observation: 


Lemma 1. The Markov chain is ergodic on f^L,Ar;ni,n 2 • More precisely, one can go from any 
V £ ^L,N-,ni,n 2 to any rj' G f^L,Ar;ni,n 2 ® chain of elementary moves where a single particle jumps 
by +ei. 


Proof of Lemma^ First let us define the height function hp^pi of r/ respective to p' m- this is 
defined on hexagonal faces / and its gradients are given by 




(3.2) 


v&C 


f^f' 


where Hri{Cf^f/) was dehned in (2.10). 

The r.h.s. of (3.2) is independent of Cf^p: this amounts to the statement that h^,T;'(/0 “ 
hri^p/(f) = 0 if f' = f, which follows from the second equality in (2.13) because both ri,ri' are in 
^L,N;ni,n 2 with the Same values of rei = Nki,n 2 = Lk 2 . The height function is therefore well- 
defined modulo an additive constant: let us fix this constant by establishing that hp^p > 0 and 
that it vanishes at some face that we call /. Given any / such that hp^p{f) > 0, construct a path 
/ = /o, /i, / 2 , • • • with the rule that the edge traversed from fi to /i+i is not occupied in rj and has 
the white vertex on the left. We see then that hp^p is non-decreasing along this path. We have: 


Lemma 2. The path f = /o, /i, / 2 , • • • cannot contain a closed loop. 


As a consequence, it must be the case that after a hnite number k of steps the path cannot be 
continued. In this case, this implies that at face fk it is possible to rotate the three dimers (all 
three edges with white vertex on the left are occupied, as in Figure]^ (&))• the vertical one moves 
by ei and the height at fk decreases by 1. Repeat the procedure until h^^p is zero everywhere. 

□ 


Proof of Lemma^ Assume by contradiction that the path contains a loop Cg^g. We know that 
for every X G f^L,N-,ni,n 2 we have 

HxiCg^g) = HgiCg^g). (3.3) 

Since all edges traversed by Cg^g have the white site on the left, Hx{Cg^g) equals minus the 
number of traversed dimers. For rj, we know that no dimer is traversed, so we conclude that 
none of the dimers traversed by Cg^g is occupied in any configuration X G IlL,Ar;ni,n 2 - This leads 
to a contradiction with our assumption ni,n 2 ,NL — ni — n 2 > 0. Suppose for instance that a 
vertical edge is contained in no A G ktL,N-,ni,n 2 - by translation invariance, this would imply rei = 0. 
Similarly, the assumption n 2 > 0 and LN — ni n 2 > 0 exclude the case where some given 
non-horizontal edge is contained in none of the conhgurations X G TlL^x-,ni,n 2 - ^ 


4. Invariance of Gibbs measures 


Our main result is: 


Theorem 1. The probability law vr := TTL,N-,mi,m 2 is stationary in time. 
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Remark 5. In [19], stationarity of t: for q = Q and all = 1 was proven, and actually it was 
shown that the infinite volume Gibbs measures obtained in the limit L —)• oo,N —>■ oo, 

mi/L —>■ pi,m 2 /N —)• p 2 , are stationary for the dynamics on the infinite hexagonal lattice. The 
proof of invariance of -k on A-l,n in the case q> is more involved than for q = 0. 

Proof of Theorem^ Call L the Markov generator, then we have to check [7r£](r^) = 0 for every 
configuration p G IlL,N-,ni,n 2 - This can be rewritten as 

'^TT{a)C{a,p) + 7r{p)C{p,p) = 0. (4.1) 

o-/’? 

Since the generator has sum zero on rows, this is equivalent to 


7r{p) 




cr 


= 0 . 


Now it is obvious from the definition of the dynamics that 

(1 - q^p){l - q^p+^) 


p 


(1 - gC'p+i) 


=: Slip). 


On the other hand, we claim (see proof below) that 
7t{p) 


^ V- (1 - q^P+^)il - q^p) ,,, 

Z^®»-(p)+i (^i_qFp+i'j • < 52 ( 1 /), 


" P 

where r(p) + 1 equals 1 if r(p) = N. Finally, we claim that 

Slip) = S 2 ip) for every p G IlL,N-,mum 2 - 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


Proof of (4.4) An update a ^ p means that a family V^fi in a has been shifted by ei. Call, as in 
rfi the highest particle in . The reverse move consists in shifting Vfi by —ei. Call 


Dehnition 


Pp the configuration obtained by p by shifting V~ by —ei. We see then that the configurations pp, 
with p running over all particles, exhaust the configurations a contributing to (4.4). Next we claim 
that 


AVp)rr ^ il - q^P+^)il - q^p) 

T^i^) ^ivp,v) ar(p}+l 


(4.6) 


If this is true, then (4.4) is proven. Label p_£+i,= ,p_£+ 2 ) ■ ■ ■ iPo ^ P the i = \Vp \ < N — 1 

particles that are shifted left in the move p —)• pp, starting from the lowest one. We have from (3.1), 

(1 - g®^’-«+i+^)(l - q^^-f+i) 


^ihp:V) 


1 - q'^P-e+1 

As for the ratio of probabilities, one sees that 

F-jpp) _ «r(po)+i {Q'iQ)Fpgiq;q) 

T^ih) ar(p_,+i) iT<l)F^^ + liT<l)c^_,^^-l 

A (9;9)Ap.+i(g;g)Dp,-i 

X M 

j=-i+i 


(4.7) 


(4.8) 


( 9 ; g)E„, (g; q) 


B„ 


iq-, q)A„^ iq-, q)D„, iq-, q)E„,-iiq; q)B„,+i 


«r(po)+i (1-g'^^-^+i) -A (1 - A^^"^")(l - 
ar(p_,+0 (1 - Ao+i) . ^ (1 _ gA )(1 - gA+1) ■ 


Ap - +1 
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Note that the F’s and C’s of intermediate particles do not appear because they are exactly zero, 


both in r] and in rjp. Also, note that Bpg = Ap^,,B. 


D 


P-e+2 ■ 


Therefore, (4.8) simplifies into 


P-f+2 


= A 


P-£+l 


and Bp_^ = D 




E 

> ^P-t+\ 


ar(po)+i (1 - - q^po+'^){l - q^po) 

Or(p_^+i) (1 - - q^^-e+i){l - q^Po+^) 


Once we multiply (4.7) by (4.9) and we recall that po = p, we obtain (4.6) and therefore (4.4). □ 


Proof of (4.5). It is sufficient to prove that Si — S 2 is independent of rj. In this case, it must 
be zero because [vrT] [rj) = 0 by conservation of the probability. Thanks to Lemma it is 
sufficient to show that Si — S 2 is unchanged when a single particle is moved by ei. The sums in 
Si, S 2 make sense also when horizontal particle positions are real numbers; given a particle p that 
can be moved by ei, we let r]{s) be the configuration where p is moved by sei,s G [0,1] and we 
prove that ds{Si — S 2 ) = 0. Note that not only the values of Ap,... ,Fp depend on s, but also 
Dp .^, Ep .^, F „3 , , Bp ^, Cpg do. We find by direct computation 


where 


dsSi 

logq 


“h ^r(p)+l*^11 ? 


Qilici ^ ^ ^r(p)+l*^21; 


1 _ 1 _ n^P 

O Sr, ■*- y -Dn + l y 

Sio = q I- C~T~1 


1 i 1 — qCp+i 

, , n„,+i 

(l_gCp+i)2 l_gCp3+l’ 

« (l_gCp3+l)2 9 1 

^ F ,l(l-(7^-3+l)(l " 

S 20 = g ^ 


gCpg+l ’ 

lp2+l 


S 2 I 


: . ... 
^ (l-g^P3+i)2 ^ i-qFp2+^' 

A„+l 1-g^^ 1 - 

- n-^^P'^ _ _ n P _ ___ 

^ 1 _ qFp+i ^ i_ qFp+i 


i'p+i ^ 2_ q^p~^^ 

-g-^r+i)(i-ggp) _ , +1 l-g^^4 

(l_^Fp+i)2 ^ l-qFp,+^ 


(4.10) 


(4.11) 


(4.12) 

(4.13) 

(4.14) 


^Pl — ^Pl — ^b-} ^Pe — ^p-J ^P6 — ^p ^pf ^PQ — -^P ^P'> ^P5 — ^P5 


Dp — Ep 


1 n^P~^^ 

Q - n^P ^ 

*^10 Q ■] _|_-| 

1 — n^p^^ 


- q^p 


n„+i 1- _ 

1 — gGp+i 1 _ qCp+i 


(4.15) 


Sii 


, „Cp+i (l-g^-)(l-g^^+^~) , . +1 ^ 

+ 9 (l_gCp+l)2 +9 


- qFpi 


1 


( 1 - 


1 ^ 


qAp—Bp-\-i 


q^^ ^^) _ E, l-g^P5+l 

qFp-\-1^2 ^ —K^-l- 


qDp Ep -\-1 


(4.16) 





























STATIONARY MEASURE OF THE q-WHITTAKER PARTICLE SYSTEM 


11 


while using Ap^ — Dp , Dp^ — Ep + Dps Dp , Fp^^ — Dp Ep , Fp^ — Cp , Ap^ — Bp + Cp , Ep^ — 

Dp Cp, Ap2 — Ap + Bp^ Dp, Ep2 — Bp, Fp2 — Ap Bp, 


5 20 = 

5 21 = 


(^1 _ ql^p+Cp+i'^^-^ _ qDp—Cp') 
(1 — gC'p+l)2 

1 _ g^P + 1 


I- qAp+ap^ 

+ qBp 


Ap-\-Bp. —Bp-\-l 


_ QAp—Bp-\-l 


1 - q^p 


- q^P- 


q 


1 — qBp+^ 1 — g^p+l 

F„+l (l -9^^+^)(l -g^P) 

— qBp^i^2 ^ 2 — qDp—Ep^i 


(4.17) 

(4.18) 


Note that, when taking the difference 5io — 5'2o, the dependence on Bp^ cancels and likewise for 
the difference 5ii — S '21 the dependence on Dp^ cancels. Finally, one checks that each of these 
differences, which a priori depends on q, Ap ,..., Ep, is actually identically zero and (4.5) is proven. 
This concludes the proof of Theorem □ 
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